Rotational stability of linear string baryon 
configuration 

V.P. Petrov, G.S. Sharov* 

Tver state university 
Tver, 170002, Sadovyj per. 35, Mathem. dep-t. 

Abstract 

For the linear baryon string model with three massive points (three 
quarks) connected sequentially by the relativistic strings the initial-boun- 
dary value problem is stated and solved in general. This problem implies 
defining a classical motion of the system on the base of given initial po- 
sition and initial velocities of string points. The given solution work 
reduces the initial-boundary value problem for the considered model to 
the system of ordinary differential equations that can be integrated nu- 
merically in general. The examples of numerical simulation of the system 
motion are considered, which, in particular, allow to answer to the old 
question about stability of a flat uniform rotation of the linear string 
configuration. 



Introduction 

String models of baryon differ from each other by geometric character of 
junction of three massive points by relativistic strings. Four variants are pos- 
sible: a) the "three-string" model or Y-configuration with three strings from 
three quarks joined in the fourth massless point EL 0|; b) the "triangle" model 
or A-configuration with pairwise connection of three quarks by three relativis- 
tic strings || (7) ; c) the quark-diquark model q-qq |§ (from the point of view of 
classical dynamics it coincides with the meson model of relativistic string with 
massive ends [Q, ||); d) the linear configuration q-q-q with quarks connected 
in series H. 

In the present work the latter model is considered. It is not investigated 
practically in comparison with others. The common opinion was that the 
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configuration q-q-q is unstable. It seemed obvious, that for the rotational 
motion of this system the centrifugal force rejects the middle quark at an 
end, and the system q-q-q transforms in quark-diquark q-qq 0. However this 
conclusion remained only speculative till now — for its proof it was necessary 
to develop any way of solution of the initial-boundary value problem for the 
given model with arbitrary initial conditions. It would allow us to investigate 
on stability, for example, the known rotational motion, in which the rectilinear 
string uniformly rotates and the middle quark is at rest at the center of rotation. 

In this work the initial-boundary value problem for classical motion of the 
linear configuration q-q-q is solved on the base of approach, applied earlier to 
the "meson" model of relativistic string with massive ends [[lO) . 

The initial-boundary value problem for the baryonic model has additional 
difficulties in comparison with the mesonic one. In particular, for an arbitrary 
motion of the q-q-q system the world surface of the string is not smooth on the 
trajectory the middle quark, moreover, there is no parametrization with the 
conformally flat induced metric and with the quark trajectories as coordinate 
lines. 

These problems and also equations of motion and boundary conditions for 
the given model are considered in Sect. 1. For solving this initial-boundary 
value problem the Fourier method (widely used in the string theory p|) is not 
applicable owing to nonlinearity of the boundary conditions and necessity of 
realization of the orthonormality conditions. Similar difficulties connected with 
the latter conditions arise in numerical solution of the problem with the help 
of difference schemes. In Sect. 2 we suggest another approach to solving the 
given problem. In section 3 with the help of numerical methods the concrete 
example of motion for the q-q-q string baryon model is investigated. Driving 
for. 



1. Equations of motion and boundary conditions 
On quark trajectories 

Let's consider an open relativistic string with the tension 7 carrying three 
pointlike masses mi, TO2, 771,3 (the masses mi and 777.3 & r e plased at the ends of 
the string). The action for this system is 

<*3(t) , , 



S[X"] = -/dr{ 7 / [{X,X') 2 -X 2 X' 2 \ 2 da + J2^X* 2 (r)\. (1) 

Here (r, a) are coordinates of a point of the string in .D-dimensional Minkowski 
space R X ' D ~ X , the speed of light c = 1, (r,cr) e fl = fli U fl 2 (Fig. 1), 



2 



(a, b) = a^bn — (pseudo)scalar product, = d T X^, X'* = d a X», X* m (t) = 
J^X m (t, cri(r)); ct^t) (i = 1, 2, 3) — inner coordinates of world lines of point- 
like masses (quarks). 
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Figure 1: Domain of integration in (1). 



We use the notations 

L(X,X') = - 7 yJ(X,X')2-X2X'Z, K{X*) 



-rn 



Varying action (1) and equating the variation 5S[X^] to zero we shall receive 
the equations of motion 



d dL 



d dL 



drdX» dodX'v 



0, 



(t,<t) e 0. 



(2) 



er = <7j(r). and boundary conditions at the quark trajectories a — CTj(r). 

To derive boundary conditions in the model "q-q-q" we are to take into 
account the discontinuities of X^, X ,fJ - on the line a = 0-2(1"). Thereby the 
term 

[VflLsX^ + ^i^SX^dTda (3) 
Ia 9t dX^ ' da K dX'^ n y ' 

in (JS'fX^ 1 ] transformed using the Green's formula equals the sum of two curvi- 
linear integrals of internal boundary values along the borders of the domains 
fii, 2 (Fig.l) and therefore — in the following boundary conditions: 



d gAi 

d^dxf 1 



dL 



dL 



dX'^ dX^ 



0, 



cr=(Ti(r) 



3 



d dA 2 
d? dX^ 



dL 



dL 



dL 



dL , 
: — &i(t) 



+ 



d dA 3 

d^axY 



dX't* dX^ 
dL dL 



a=cr 2 (T)+0 

= o, 

cr=(T2(r)— 



(4) 



dX'v dXf+ 



0, 



a=a 3 ( T ) 



The equations of motion of the string and the boundary conditions have the 
simplest form if with the help of nondegenerate reparametrization r = r(f , a), 
a = a(f,a) the induced metric on the world surface of the string is made 
continuous and conformally-flat H,i.e., satisfies the orthonormality conditions. 



X 2 + X' 2 =0 7 (X,X')=Q. 
Under conditions (5) the equations of motion (2) become linear 

X» - X"' 1 = 
and the boundary conditions (4) take the simplest form 



(5) 
(6) 



m 1 -U 1 ^T)- 1 (X^ + a[(r)X^ 



0. 



(7=<Ti(t) 



m 2 —U2(T) +1 (X'» + a' 2 (T)X») 



7(X'^+4(r)^ 



ct=<t 2 (t)-0 

= o, 

a— <T2 (t)+0 



(7) 



m 3 -U^r) +1 (X^ + a' 3 (T)X^ 



0. 



0-=o- 3 (t) 



Here the notations 



U?(t) 



VWi^) Jfr . (1 - ^) 



1, 2,3 



(8) 



for unit i? 1 ^- D ~ 1 -velocity vector of i-th quark was taken. 

From the physical point of view equations (4) are the 2-nd Newtonian law 
for the material points m,. 

Eqs. (5) and (6) are invariant with respect to reparametrizations r ± a = 
f± (f ± a) [9 . Choosing these two arbitrary functions f± one can fix two (of 
three) functions crj(r), for example, in the form 



cti(t) = 0, Cr 3 (T)=7T. 



(9) 
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The first equation (9) may be obtained at the first step by using the above 
reparametrization with required /+ and f-(rj) = r), the second equation (9) 
- at the second step with /+ = /_. This procedure doesn't permit to fix 
ct 2 = const for all r in general. 

For the initial-boundary value problem the function <J2{j) should be is 
calculated on the base of initial data. 

The function (r, a) is continuous in f2 but on the line cr 2 ( T ) its derivatives 
can have discontinuities of 1-st type (except for tangential and -§^Vt)- 

2. Initial-boundary value problem 

The initial-boundary value problem for a relativistic string with the action 
(1) we shall formulate in general as follows: " to find a solution of the equation 
(6) X^(t, a), fj, = 0, 1, . . . , D — 1, smooth in the domain 

fi = (r,cr) : r > to(o-); <7 X (t) < a < cr 3 (r). (10) 

continuously differentiable on dfl satisfying to the orthonormality conditions 
(5) and boundary conditions (7) and two given initial conditions: an initial 
position of the string in Minkowski space and an initial velocity of a point of 
the string" . 

An initial position of the string can be given as a spacelike curve in Minkowski 
space 

^=p"(A), A G [Ai, A3], (11) 

satisfying the conditions [//(A)] 2 < 0. Initial velocities on the initial curve can 
be given by timelike R 1 ' ' 1 - vector v^(X), A G [Ai,A 3 ]. 

To solve the problem we set the initial curve (11) on the world surface of 
the string parametrically 

t = t(A), a = a(X), Ag[Ai,A 3 ] 

(in expression (10) the same line is given as r = To (a)). 

In these notations the initial conditions (initial position of a string and the 
initial velocities) can be given in general jl(J |ll| 

X»(t(X),<t(X)) =p"(A), Ai<A<A 3 , (12) 
a(X) X"(t(A),«t(A)) + /3(A) X^(r(X),a(X)) = p»(X), A x < A < A3, (13) 

where a(A), /3(A) are arbitrary functions satisfying the inequality a(X) > |/3(A)| 
and (if the relations (9) are tied) the condition on at the ends /3(Ai) = /3(A 3 ) = 
0. 

Functions t(A), <t(A), ot(X), /3(A) are related by formulas flicj 
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,,,, aA + pP , /3A + aP 

* W = r^—> A*) = 3 — > ( 14 ) 



where P(A) = (v(X), p'(X)) , A (A) = y (u, p'} 2 — w 2 p' 2 . It is convenient to 
choose the constants of integration in (14) so that r(Ai) = <r(Ai) = (Fig. 2). 
Realization of the condition ct(A 3 ) = ir is possible by multiplying a(A) and 
/3(A) by a constant. 




The line (72 (t) devides domains f2 into two domains f2i and f^- 
consider the common solutions of the equation (6) in these domains 

X»(t, a) = \ [*? + (t + a) + *?_(t - a)] , (r, a) e fli 
X»(t, a) = \ [^ + (t + a) + *£_(t - a)] , (r, a) e fi 2 . 



We can find the functions "J" 



'(A) ± (7(A)) 



1, 2 using the formulas |10 

(A =f P)v» ± v 2 p'» 
A(a(A)±/3(A)) 



If (as mentioned above) we choose 

t-(Ai) = o-(Ai) = 0, 



<7(A 3 ) = 7T, 



Let's 



(15) 



(16) 
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that the formulas (16) will determine ^y_(^) in the following segments: 

£ G [0, r(A 2 ) + a(A 2 )] , £ G [r(A 2 ) + a(A 2 ), r(A 3 ) + tt] , 

*?_(0, ^G [r(A 2 )-a(A 2 ),0], *£_(£), £e [r(A 3 ) - tt,t(A 2 ) - a(A 2 )] , 

(17) 

that allows us to find the solution of the problem as (15) in the zones A = 
{(r,a) : r(A(a)) < r < <j;t + a < r(A 2 ) +a(A 2 )} B = {(r,a) : t(A(ct)) < 
t < a + t(A 2 ) — er(A 2 ); t + a < t(A 3 ) + 7r} (Fig. 2). In these zones the solution 
depends only on initial data without influence of the boundaries. The constant 
of the integration in Eq. (16) is determined with the help of the initial condition 
(12). 

In others zones of the domains fii and Vt 2 we shall Obtain the solution by 
prolongating the functions "J/y. on all positive semiaxis with the help of the 
boundary conditions. 

For the given initial-boundary value problem the boundary conditions look 
like (7). By choosing the inner equations of the string ends as (9) and by 
substituting in Eqs. (7) common solutions (15) we transform this system (it is 
supposed that the masses mi, m 2 , m 3 are final and nonzero): 



dr 2mi 



[K + (r)-K-(r)]=0, (18) 
CI ^a' 2 )[^ + (r + a 2 )-^ + (r + a 2 )} + 

= 0, 



dr 2m 2 

+ (1 - <7 2 ) [*£_(t - a 2 ) - *y_(T - a 2 )] 



(19) 



^ + -2_[^ +(T + 7r) -^_ (T - 7r)]= 0. 



(20) 



The expressions for the unit velocity vectors (8) after substitutions (9), (16) 
may be rewritten as 

mT)= *'iW + *T-(t) _ ^ + (r + 7r) + *^_(r-7r) 



(1 + a> 2 ) *£(t + a 2 ) + (1 - a 2 ) *£_(t a 2 ) 



= (l + 4) + a 2 ) + (l - <y' 2 ) *t(r - Q2) (21) 

v /2(l-a 2 2 )(* 2+ (r + ( r 2 ),vl/ 2 _(r-a 2 )) 

In last expression the continuity of the function X^{r,a) and its tangential 
derivatives on the line a = cr 2 (r) was taken into account. By differentiating on 
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t the identity 



(7 — (72 (t)+0 



(after substition (15)) we obtain 



(l + ^)^ + (+) + (l-^)^_(-) = (l + ^)^ + (+) + (l-^)^_(-). (22) 



Here the notations (+) = (r + CT2(t)), (— ) = (r — 02(1")). 

We transform the systems of the ordinary differential equations (18) - (21) to 
the normal form by the method used in Ref ]l0| . Let's consider transformation 
of the boundary condition (19) with the help given method. 

Using the scalar product of the Eq. (21) by (1 + a' 2 ) *■+(+) and (1 - 
<t' 2 ) ^'IL{— ) and taking into account the isotropy of the vectors 9^. 



^ y J{l-a^)(% + (+), *!_(-)> = (1 + <j' 2 ){U 2 (t), = 
= (l-a' 2 )(U 2 (T), %_(-)) = ^{l-a'i)(% + {+),%_{-))= (23) 
= (l + a' 2 )(U 2 (r),% + (+)) = (1-*' 2 )(U 2 (t),%_(-)). 



(1 + 4) *&(+) + (1 - 4) *£H = 2(1 ± ^)<C/ 2 (r), *; ± (±))t^(r), 



(right for all 8 variants of a choice of indexes i, j and signs ±) with the help 
of which we shall transform the summands in Eq. (19) as follows: 



(i + a' 2 ) - (i - ^2 ) *£.(-) = 2(1 - a> 2 ) [(u 2 , *U-K *T-H] , 

(1 - (-) - (1 + ^) * 2 M +(+) = 2(1 + <r£) [<E/ 2 , # 2+ (+)>t/£ - . 



= JL [jm _ ^ (r) [ (1 + ^) (+) + (1 _ 4) , (24) 



Using the same receptions we shall transform the equations (18) and (20) 
to the form 




The equalities (21) and (23) result in the relations 





[«-^W^(r)]<(r), 



(25) 



g?t mi 



dUjij) _ 7 



[5Z-Ug(T)U Sv (T)]%»_{ 



(26) 



dr TO3 
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The systems (24) - (26) need the initial conditions 



C/f(r(A 4 )) -^(AO/v^W. i = l,2,3. (27) 

Integrating systems (24) - (26) with the initial conditions (27) we can de- 
termine for r > r(Aj) unknown vector functions Uj* (r) with the help of the 
functions "J/-^. known from the initial data on the segments (17). This proce- 
dure is limited on r to ordinates of points in which the trajectories a = Uj(r), 
i = 1, 2, 3 cross the characteristic lines r ± a — const (Fig. 2). However we can 
continue this procedure for all r if for every value of r we determine the func- 
tions ^ outside segments (17) with the help of Eqs. (18) -(20), (24)- (26). 

For this purpose we rewrite Eqs. (18) and (20) as 

*;i(r) = *'4(r) - 2m l7 - 1 ^(r), (28) 
*&(t + tt) = *£_(r - tt) - 2m 37 - 1 ^(r), (29) 

Eq. (28) corresponds to Eq. (25) and Eq. (29) — to Eq. (26). From Eqs. (19) 
and (22) we deduce two similar relations 

*£-(t - <T2(t)) = *?L(t - <7 2 (r)) - ma [ 7 (1 - ^)]"V^(r). 

For solving the system (24) we are to determine the function <j 2 (t) for 
t > t(X 2 ). This function can be obtained by taking two of the equalities (23) 

(1 - a' 2 ){U 2 (r), *;_(-)) = (1 + a' 2 )(U 2 (r), *' 2+ (+)) 

and expressing a 2 as follows: 

da 2 (r) (U 2 (r),[%_(r-a 2 )-% + (T + a 2 )}) 
dr {U 2 {T),[V l _{T-a 2 ) + % + (T + a 2 )]y 1 > 

The equations and systems (24) -(31) allow us to continue the functions 
determined from the initial conditions with the help of formulas (16) on 
the segments (17). The constants of integration in are determined from 
the initial condition (12). The descibed above algorithm of calculating 
(with taking into account Eq. (15)) is the algorithm of solving the considered 
initial-boundary value problem with arbitrary initial conditions p M (A), v^(X). 

3. Numerical solution of initial-boundary value problem 

To the present time only two elementary motions of string configurations 
q-q-q are known. They are analytical solutions of equations (6) and satisfy 
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conditions (5), (7). One of these motions is spatially one-dimensional com- 
pression of a symmetric (mi —1113, vn-i is at the center) rectilinear string that 
in some initial instant was at rest S. The second motion is a flat uniform 
rotation of the rectilinear string with the middle quark at rest at a center of 
rotation (in the reference frame of this center) || ^, || . If the choice of r, a 
satisfies the conditions (5) and (9) then this solution can be presented as 

X^ — W |ar; sin(acr — <Pq) ■ cos ar; sin(aer — 4>o) • sin err}. (32) 

Three coordinates of the vector X^ are X° — t, X 1 = x, X 2 = y. The 
trajectory of the middle quark on a plane r, a has the form ct 2 (t) = <p a /a = 
const and the trajectories of the others quarks have the form (9). 

The vector function (32) is the solution of the equations (6), satisfies the 
conditions (5) and under the relations 

7 sin 0o 7 sin(a7r-0o) , , > ^ n 

cos(acTi - <po) > (33) 



Lumi cos 2 0o ' urn?, cos 2 (a7r — 0o) 

it satisfies the boundary conditions (7). 

The solution (32) describes a rotation of the rectilinear string with an an- 
gular frequency to. The quarks move at constant speeds 

Vi = |sin(a<7i - O )|, i = 1,2,3 (34) 

along the circles with the radii Ri = Vi/uj. The relations (33) may be rewritten 

as 

= 1 2 = r> r-i 2T' * = 1, 3. (35) 

mi 1-vf Ri(l-vf) 

The expression (32) can be obtained as the solution of the initial-boundary 
value problem (6), (7), (12), (13) with the initial conditions 

p^(A) = {0;-i? 1 + ( J R 1 + J R 3 )A;0}, _ Rl 

^(A) = {1;0;-, 1 + ( U1 +, 3 )A}, AG[ °' 1J ' ^ " Jfe+ife- (36) 

If the parameters Uj, Ri satisfy Eqs. (35) the solution of the initial-boundary 
value problem has the form (32) in all domain and can be obtained with the 
help of algorithm described in the previous section. 

To solve the problem of stability of the solution (32) we are to consider an 
initial-boundary value problem with the initial conditions p^(X) w M (A) which 
slightly differ from conditions (36). In particular, it is possible to put A2 ^ 
i?i/(-Ri+i?3) = Vi/(vi+vs), i.e. to place the middle quark not at the rotational 
center. 
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Figure 3: Motion of the system with initial conditions (36), (37). At = 0.15. 

For arbitrary initial conditions the solution of the problem becomes com- 
plicated. In particular, system (24) -(26), (31) are solvable analytically only 
for exclusive cases. Therefore the use of numerical method is necessary. It 
makes possible to solve the initial-boundary value problem with any initial 
and boundary conditions. 

The used below procedure of solving includes three stages: 

1) realization of initial data, i.e., the determinition the functions *E^j_ in the 
segments (17) by the formulas (16) [with the help the initial conditions //'(A) 

«"(*)]• 

2) continuating the functions ^i±(0 (with integrating) with the help of 
boundary conditions (24) -(31); 

3) searching of level lines of the world surface t — X°(t,(t) = const and 
drawing of their projections on the plane Oxy (x = X 1 , y = X 2 ). 

The last item is required for the best representation of the string motion. 
We show " photographs" of its sequential positions made at regular intervals 
At. 

The results of calculations for the q-q-q system with mi = 1, mi = 3, 
7713 = 2,7=1 obtained with use of the MATLAB package are shown in Fig. 3. 
The initial conditions for the given motion have form close to (36), where 

ui = 2, wi ~ 0.4142, v 3 ~ 0.4142, A 2 = 0.615 ^ v 1 /(vl + v3) ~ 0.637. 

(37) 

The values v\ satisfy conditions (35) however the middle quark is located 
not at the center of rotation. 

In Figs. 3a- d on the plane Oxy the positions of the system are shown 
sequentially with the interval At = 0.15 in time. The initial position of the 
string (36) in Fig. 3 is marked by number 1 and next (with growth of t) "pho- 
tographs" of the system are marked by numbers 2,4,..., near the third quark 
position. For convenience of perception the first pointlike quark q\ with mass 
mi = 1 at the end of the string is depicted as a green square, the third quark 
g3 with mass 1113 = 2 — as a red square, and positions of the middle quark g 2 
are marked by blue circles. 

The last position of the system number 11 in Fig. 3 is the first one in Fig. 3b. 
The system q-q-q continues rotating as a whole in the counter-clockwise direc- 
tion in Figs. 3b- d. At the beginning (Fig. 3) the motion of the system as a 
whole is close to the uniform rotation (32). Later the middle quark qi under 
the centrifugal force is displaced to the quark q\ but there is no junction of qi 
with q\ and transformation of the system in the quark-diquark one p|. The 
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middle quark plays a role of rotational center for the segment of the string 
92-91 (positions 26, 27 in Fig. 3c) and then it returns to the center of a system 
(Fig. 3c, 3d). The middle quark does not reach this center and begins again 
centrifugal moving away with recurring of all cycle. 

As we see, the motion is quasiperiodic. It's qualitative features are the same 
for rather wide variation of parameters of the initial-boundary value problems. 
In particular, decreasing of the initial deviation of 92 from the center of masses 
reduces in retarding of the initial stage of motion. 

If the initial conditions differ from (36) by a small nonsymmetrical cor- 
rection to the velocity v^(\) (for A 2 = R ^ Ra ), the motion of the system is 
approximately the same as in Fig. 3. The increase of the mass mi of the middle 
quark q2 reduces in Increase of AR — minimal length of the segment q\-q2- 

Energy or velocity of rotation of the system under the relations (35) weakly 
influences to the character of motion. In the non-relativistic limit Vi -C 1 and 
in relativistic case (Fig. 3) the period of cyclical motion of the middle quark 
between the center and the end is approximately equal to the rotational time 
of the string. 



Conclusion 

In the present work the method of solution of the initial-boundary value 
problem for classical motion of the linear string baryon model q-q-q is sug- 
gested. This approach allows us to solve the stability problem for the known 
solution (32) — flat uniform rotation of the rectilinear string with the quark 
92 at rest. This solution is unstable in Lyapunov's sense — the as much as 
small perturbations result in complicated motion (Fig. 3) with quasi-periodical 
varying of the distance between the nearest two quarks. However the minimal 
value of the mentioned distance AR does not equal zero, in other words, the 
system q-q-q is not transformed in quark-diquark (q-qq) one, as was supposed 
speculatively j|. 

The analysis of this problem does not permit us to exclude the linear con- 
figuration q-q-q from a set of string baryon configurations (q-qq, Y, A) for 
describing baryonic orbital excitations on the Regge trajectories. At that it is 
possible to use expressions for the total energy and the angular moment of the 
system q-q-q, obtained in Ref. [0 for the motion (32). 



3 r 



3 
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Here AE,s is the spin-orbital correction to the energy, s, are projections of 
quark spins, V2 = 0. These expressions are applicable and for more compli- 
cated motions, obtained with the help of weak perturbations of the conditions 
(36) (see Fig. 3) because of conservation of E and J. Hence, for perturbed 
quasiperiodic motions the graph of dependence J = J(E 2 ) is close to linear, 
that allows us to describe baryonic Regge trajectories with the help of the 
model q-q-q 
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